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Abstract. As it is well known the topology of space is not totally determined by 
Einstein's equations. It is considered a massless scalar quantum field in a static 
Euclidean space of dimension 3. The expectation value for the energy density in all 
compact orientable Euclidean 3-spaces are obtained in this work as a finite summation 
of Epstein type zeta functions. The Casimir energy density for these particular 
manifolds is independent of the type of coupling with curvature. A numerical plot 
of the result inside each Dirichlct region is obtained. 



PACS numbers: 04.62. +v; 03.70.+k 



Casimir Effect in E 3 closed spaces 



2 



1. Introduction 

The Universe today is very close to homogeneous and isotropic. Then the most 
important 3-Manifolds in connection with cosmology are S 3 , H 3 and E 3 . It is well 
known that R x S 3 , R x H 3 and R x E 3 are solutions of Einstein's equations (EQ). It 
is irrelevant for the EQ, if the spatial sections have a non trivial topology modelled on 
them, or not. 

There are many effects that arise if space is compact. The first astrophysical limits 
on the topology of the universe were obtained for a 3-torus T 3 . Accordance with the 
homogeneity of the cosmic microwave background radiation CMBR puts a lower limit 
on the size of the fundamental cell, about 3000 Mpc, which is a cube in the cases 
of p] and [2]. Later on, it was shown that this result is very sensitive to the type 
of the compactifications of the spatial sections. For a universe with spatial sections 
T 2 x R, the fundamental cell's size is about 1/10 of the horizon, and is compatible 
with the homogeneity of the CMBR [3]. In compact universes, the pair separation 
histogram would present spikes for characteristic distances. At first it was thought that 
this technique, known as the crystallographic method, was able to totally determine the 
topology of the universe pE]. It turned out that the crystallographic method only applied 
when the group that defines the manifold, contained at least a Clifford translation, i.e. a 
translation which moves all the points by the same distance [5] and [6] . Generalisations 
of the crystallographic method were proposed, for example in [TJ. 

Also in compact universes the light front of the CMBR interacts with itself 
producing circles in its sky pattern [8]. 

Recent results have called our attention to the possibility that methods based on 
multiple images will prove not to be efficient |9j. The reasoning is that, according to 
observations, the curvature is very small, so the fundamental regions are so big that 
there has not been time enough for the formation of ghost images. The result is that for 
low curvature universes such as ours, only compact universes with the smallest volumes 
could be detected by pattern repetitions. 

The knowledge of the spectrum of the Laplacian for compact spaces have some 
motivations. One of them is the decomposition of the CMBR in right set of modes, if 
space is non trivial. It has been studied numerically for compact hyperbolic space by 
[TOf [TT]. And more recently [12] for compact spherical space. 

In 1948, Casimir calculated an attraction force between parallel non charged plates 
in vacuum [13J. This force is universal, in the sense that it does not depend on the 
physical properties of the plates, or the electric charge. Its origin is of topological 
nature, connected to the boundary conditions. There are interesting articles, in which 
the imposition of boundary conditions is addressed in a more rigorous and realistic 
approach, for instance [14"] . 

The Casimir effect also occurs in compact spaces. For the present era where the 
size of the Universe is ~ 3000 Mpc the Casimir effect is irrelevant. Anyway, soon after 
the cosmological singularity the co-moving distances were on the order of the Planck 
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scale. It is general believed that the very primordial Universe was dominated by strict 
quantum effects including Casimir effect. Although today quantum effects are commonly 
disregarded, they must have been fundamentally important in the primordial Universe. 
We shall use the massless scalar field toy model. 

The Casimir effect can be obtained by analytic continuation of the zeta function. 
This procedure is known as the zeta function regularisation, see for example [15], |16j . 
Anyway, in a more realistic context, the ability of the zeta function method has been 
questioned [17] . For a calculation of the Casimir energy in the spherical space, see [18J. 

In this work we expose a different technique which was previously used by us for 
the hyperbolic space [T9l [201 [2T] . 

Previous calculations for the Casimir effect for the Euclidean space were done, for 
instance by, DeWitt, Hart, Isham and others [22]. Particularly interesting is Dowker's 
article in which the result is also written in terms of Epstein's zeta functions. A 
topological classification of possible "gauge fields" over compact base manifolds, is also 
given in [22]. These and other very interesting topological questions are not going to be 
addressed in this present work. 

The crystallographic groups are discrete subgroups of the full isometry group of the 
Euclidean space. Among all crystallographic groups, only a few of them are of interest 
in connection to manifolds. These groups must be torsion-free, for the Euclidean and 
Hyperbolic spaces. In the celebrated list of problems proposed by Hilbert, the 18th was 
answered affirmatively by L. Bieberbach in 1910-1912 [23], [24] , see also [25]. Bieberbach 
showed that there are only a finite number of crystallographic groups for an Euclidean 
space of fixed dimension. For example, when the dimension of the Euclidean space is 3 
there are only 10 compact manifolds, of which 6 are orientable and 4 non orientable. 

In this present work we shall investigate the 6 compact orientable manifolds and 
obtain the Casimir energy for each of them. Due to the more simple structure connected 
to these topological spaces compared to the hyperbolic or spherical case, an analytical 
result is obtained in terms of Epstein zeta functions. It is explicitly obtained that 
our result is independent of the type of coupling with curvature £ for these particular 
manifolds, which means that the calculated Casimir energy density is the same even in 
an inflationary de Sitter regime. 

The paper is organised as follows. Section [2] is included for completeness. In this 
section, the expectation value for the energy momentum tensor, for a massless scalar 
field in closed Minkowski spaces is obtained. Section [3] presents a brief description of 
the compact flat manifolds. The Casimir energy is obtained in each subsection, 13.1113.61 
and a numerical plot is given. Natural units are used, c = G = h = 1, and the metric 
V = [-1,1,1,1] 

2. Casimir effect in closed spaces 

The point splitting method was constructed to obtain the renormalised (finite) 
expectation values for the quantum mechanical operators. It is based on the Schwinger 
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formalism [26] , and was developed in the context of curved space by DeWitt [27J . Further 
details are contained in the articles by Christensen [28], [29] . For a review, see [31] . 
Metric variations in the scalar action 

s = -\j ^9^A P + + m , 

give the classical energy-momentum tensor 

T, u = (1 - 20<M,, - Q " 2£j 4>A P 9,u - 2£# ;M „ 

+ 2f^D0 + £G^ 2 - V</> 2 , (1) 

where is the Einstein tensor. As expected, for conformal coupling, £ = 1/6 with 
massless fields, it can be verified that the trace of the above tensor is identically zero 
if m = 0. In the following we shall restrict to spatial Euclidean sections and massless 
fields. 

The renormalised energy-momentum tensor involves field products at the same 
space-time point. Thus the idea is to calculate the averaged products at separate points, 
x and x', taking the limit x' — > x in the end. 

(0\T flv (x)\0) = -\imT(x,x') IMU (2) 

Z x — >x 

where for the particular manifolds addressed in this article, 

T(x,x') ltu = (3) 

(1 - 2£)({d ll d l/ + d^d u ) - 2^d al + d a d p )} G$(x, x') (4) 

where is the Hadamard function, which is the expectation value of the anti- 
commutator of <f)(x) and <fi(x') on the manifold A4, see below. 

The causal Green function or Feynman propagator, for the infinite covering space 
R 3 , is obtained as 

G{x,x') = i(0\T(f)(x)(f)(x')\0} , 

where T is the time-ordering operator. Taking its real and imaginary parts, 

G(x, x') = G s (x, x') + x') , (5) 

we get for the Hadamard function 

(£,£') = (0\{<f){x),<f){x')}\0) =21mG{x,x') . 

For a massless scalar field 

GV>(x x') = — - = — - 

K x,xj 2n 2 _( t _ f y + r 2 27i 2 {x-x') 2 ' 

Green functions, as any other function defined in the spatially compact space-time 
R x Ai, must have the same periodicities of the manifold M. itself. One way of imposing 
this periodicity is by the direct summation of all the elements of the group T which 
defines the manifold 



Casimir Effect in E 3 closed spaces 



5 



The above expression is named the Poincare series, and when it converges, it defines 
functions on the manifold Ai, in this particular case, the Hadamard function. 

There are recent articles were the physical context connected to particular boundary 
conditions, Dirichlet, Neumann and periodic is addressed [32]. Edery, develops a very 
interesting cut-off technique and obtain analytical formulas as sums over gamma and 
Riemann zeta functions plus an exact multidimensional remainder expressed as sums 
over modified Bessel functions. 

Anyway if the spatial sections are compact, the point x and the point 7 a; where 
7 G T, x = 7X for all the elements of the fundamental group V. That is, the point x 
and jx are the same and identical point. This means that the only possibility is that 
the function ([6]) must be periodic in some sense. Mathematically, the functions defined 
on a closed manifold M. are called automorfic, see for example [33 J. 

In [32], for instance, the summations are over the spectrum of the Laplace operator. 
It is well known that summations over the spectrum are equivalent to summation over 
the closed geodesies, also known as the method of images, [33]. This result lies at the 
hart of the Selberg formalism for the calculations of functional traces and is valid at 
least for manifolds, Lie groups included [31], see also [T6] 

Wald's axioms for the expectation value of the energy momentum tensor are 
satisfied [35]. Albeit this does not prove that the chosen state is the correct one, it 
does indeed indicates a criterion for a candidate physical state. 

The Casimir energy density is independent of £ and is given by the 00 component 
of (jl]) which is the projection to static observers (0\T^ (x) \0)u^u u with 4 velocity 



where the / means that the direct path with 7 = 1 is to be excluded. It can be shown 
that the exclusion of the direct path corresponds to a renormalisation of the cosmological 
constant [20] . 

3. The closed oriented Euclidean spaces 

Differently from the hyperbolic or spherical spaces, only a finite number of topologies 
can be modelled in the Euclidean space. A representation for an arbitrary element for 
the rigid motion subgroup of the full isometry group can be given 



u = (1,0,0,0) 




(7) 



/ 




R 



(8) 



7 








1 / 
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where R G 50(3) is an orthogonal rotation matrix and a is the translation. The action 
of the group on any point x = {xi, X2, £3} of E 3 is linear and multiplicative 



/ 



7X 



R 







a 2 
1 



\ 



/ Xi 

V 1 / 



The discrete subgroups of the isometry group define the manifolds modelled in the 
space of interest. A group T is said to be torsion- free if for each element 7 G T, 7™ 7^ 1 for 
any integer n. In the Euclidean case, the torsion-free elements must be "screw motions" , 
i.e. a combined rotation and translation along the rotation axis. The holonomy group 
is the rotational part of T, \1/ = T/T, where T is the translational part of the group. 
Two oriented manifolds are diffeomorfic if and only if the discrete torsion-free groups 
have isomorfic holonomy groups, Of course, ^ C 50(3) is a finite group. According 
to Lemma 3.5.4 of Wolf's book [23], for the E 3 case, the only possibilities for \1/ are 
v|/ = Zi, Z2, Z4, Z3, Z6, Z2 ® Z2, which give rise to to the 6 closed oriented manifolds [23J. 
These manifolds and the Casimir energy are given in the subsections I3.1H3.6I 

The very important Theorem 8.2.5 of Ratcliffe's book [21] is reproduced in the 
following 



Every compact, n- dimensional, Euclidean space-form is finitely covered by an Eu- 
clidean n-torus. 



This theorem is in fact Bieberbach's first theorem, which states that T is a normal 
subgroup of finite index in V. This finite index is given by the order of ^ and corresponds 
to the number of coverings of the space in question by the Euclidean n-torus. The above 
theorem is used to rewrite the summations in the expectation values over all the elements 
of the group as a finite sum of Epstein type zeta functions. The Epstein zeta function 
is already well known to give the appropriate result for the Casimir energy for the n— 
torus, see for example [30] [IE]. In the original Epstein's notation let p be a positive 
integer and 

9 = (9i,-,9 P ), ft^R, 9 G R p 
h = (hi, ...,h p ), hi G R, h G R p 
m = (mi, ...,m p ), m; 6 Z, rh G Z p . 

The scalar product is defined as 

v 

(g,h) = ^9uK, 

and c^j, an invertible p x p matrix associated with a quadratic form as follows 

p 

(p(x) c^ u x ^x u . 
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Figure 1. Dirichlet Region 



With the above definitions Epstein's zeta function is written as 



9i 
hi 



9p 
h„ 



J2 '{v{m + g)]-ie 27ri ^' fl \ 



(9) 



mi,...,m p =— oo 

where the / means that mi = 0, m 2 = 0, ...m p = is to be excluded. Generalisations 
for non homogeneous Epstein zeta functions and reflection formulae, can be found in 

EE], El. 



3.1. * = Zi 

The volume of this manifold is V = ab 2 c 3 . The Dirichlet region is the parallelepiped 
given in Figure [1] According to the representation ([51) , the matrices 



7i 



73 



/ 


1 








a i 




/ 


1 








h 


\ 







i 















1 





b 2 












1 





, 72 = 










1 







V 











l ) 




V 











1 


/ 


/ 


i 



























1 





c 2 






















1 


c 3 














V 











1 , 















(10) 



with the respective inverses are also the generators of the translational part of the group 
r = T = {71,72,73}- The identifications are as follows 



r = [x,y,z) 
and according to ([7]) 



(x + ma + nb\ + kci, y + nb 2 + kc 2 , z + kc^ 



^ +00 

o% x = — 2 ' { K ma + nbl + kci ' 

m,n,k=—oo 

+ (nb 2 + kc 2 ) 2 + (kc 3 ) 2 }- 2 } 
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(4) 



where according to ([9]) the matrix 

^ a 2 abi 
c^u = abi b\ + b\ 
y aci b\C\ + b 2 c 2 

defines the positive quadratic form ip{m 
2c 23 nk. 



ac\ 
b\Ci + b 2 c 2 

cl + 4 + 4 ) 



(ii) 



Cum 2 + c 22 n 2 + c^k 2 + 2ci 2 mn + Ic^mk + 



3.2. ^ = Z 9 



The volume of this manifold is V = ab 2 c Figure [2 Now the matrices and their inverses 



/ 



7i 



1 




/ -1 



7s 



V 



o 
o 








1 




-1 






a \ 



1 



72 



/ 1 








1 



/ 



o\ 





(12) 



generate the group Y. 

The generators of the translational part T, now, are given by T = T/Z 2 = 
{71572,73}, together with the inverses. The fundamental domain, FD, for TfZ 2 has 
volume V = 1ab 2 c and is given by two copies of the Dirichlet region of Figure [2j 
FD = {l.DR, l3 .DR}. 
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-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 i 



X 

Figure 3. The Casimir energy density in natural units, eq. fTB")) inside the DR. in 
Figure [21 The z coordinate is chosen to be z = 0.6 



The identifications in this manifold are provided by 

r = (x, y, z) — > (x + ma + n&i, y + nb 2 , z + 2kc) 
— > (—x + ma + nb\, —y + n& 2 , z + (2k + l)c) 
Using (17j) , the Casimir energy is given by 

+oo 

^ 2 = -- ^ '{[(ma + n6 1 ) 2 + (^ 2 ) 2 + (2A ; c)T 2 

m,n,k=— oo 

+ [(-2a; + ma + n 0l ) 2 + (-2y + nb 2 f + ((2A: + l)c) 2 ] _2 } 






ffx 9y 9z 



(4) v + ^ 








(4), 



where 



2x 2ybi 

a ab 2 ' 

a 2 a&i 

afei 6 2 + b\ 



2y 



1/2 



\ 





V 











4c 2 



/ 



(13) 

(14) 
(15) 



The values of a = 1.2, b\ = 0.7, b 2 = 1.4 and c = 1.3 are chosen. The truncation in (Q 
is done such that the summation is over a geodesic sphere of radius R ~ 100 of integer 
numbers. The corresponding truncation error on the energy density should be less than 
1%. The result is shown in Figure [31 and is independent of the coordinate z. 



3.3. * = Z 4 

The volume of this manifold is V — a 2 b Figure HI The matrices 
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y 




Figure 4. Dirichlet Region. The opposite faces in the z direction are identified with 
a rotation of 7r/2 



71 



73 



/ 1 





/o 
1 





/ 









a \ 



1 

o\ 



b 



72 





1 






\ 



a 

1 



(16) 



together with their inverses are the generators of 17 The translational part is generated 
by T = T/Z 4 = {71,72,73}, together with the inverses. The fundamental region, 
FR for T/Z 4 has volume V = 4a 2 b and is given by 4 copies of the DR in Figure HJ 
ER = {DR,'~f 3 .DR,'~fl.DR,'~f^.DR}. The identifications and the direct summation are 



given in the Appendix A.l The result for the Casimir energy according to ([7]) and 
is the following 



7T- 








(4), 







9x 9 y 9 



(4), 



+Z 





9l 9y 9z 



(4), 



Z 





9l 9y 9l 



(4), 



(17) 



where 



x+y x-y 
a ' a ' 
-2x -2y 



2/4 



a a 
y — x x + y 



a a 

( a 2 \ 

a 2 

16b 2 



3/4 
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The values of a = 1.2, b = 1.4 are chosen. The truncation in (Q is done such that 
the summation is over a geodesic sphere of radius R ~ 100 of integer numbers. The 
corresponding truncation error on the energy density should be less than 1%. The result 
is shown in Figure El and is independent of the coordinate z 

3.4. * = Z 3 

The volume of this manifold is V — ^/3a 2 b. The Dirichlet domain is an hexagonal prism 
of side a/v^3 Figure [6] The face pairing matrices 

/ 1 -a/2 \ 

1 V3a/2 

1 

v 1 ) 



7i 



/ 1 a \ 

10 

10 

V 1 / 



72 
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7s 



/ 1 






74 



/ -1/2 


V 



a/2 \ 
V3a/2 

1 

-V3/2 
-1/2 









b 



(19) 



together with the inverses generate the group T. 

The generators of the translational part are given by T = T/Z 3 = {71, 73, 74}, with 
the inverses. The fundamental region, FR for T/Z 3 has volume V = 3\^3a 2 b and is given 
by 3 copies of the DR in Figure EJ FR = {DR^^.DR^f.DR}. The identifications are 
given in the Appendix A. 2 The Casimir energy, according to ([7]) and (Q is 



7T Z 








(4)*, + Z 





gl g] gl 



(4), 



+z 



000 

gl g 2 y gl 



(4), 



(20) 



where 



-2x x — y\/3 



1/3 



x 



yV3 x + yV3 



2/3 



( ^ 

a 2 /2 

V 







a 2 / 2 
a 2 







9b 2 



(21) 



The values of a = 1.2, b = 1.4 are chosen. The truncation in (Q is done such that 
the summation is over a geodesic sphere of radius R ~ 100 of integer numbers. The 
corresponding truncation error on the energy density should be less than 1%. The result 
is shown in Figure d and is independent of the z coordinate. 



3.5. ^ = Z 6 

The Dirichlet domain is a prism with hexagonal base of side a / -\/3, Figure [8] The volume 
of this manifold is V = \/3a 2 6. The face pairing matrices are 



7i = 



/ 


1 








a 


\ 




/ 


1 








-a/2 


\ 







1 

















1 





V3a/2 












1 







, 72 = 










1 







\ 











1 


/ 




V 











1 


/ 
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Figure 7. The Casimir energy density in natural units, given in (|20p inside the DR. 
in Figure [5] with the coordinate z = 0.6. The region is a disc of radius r = 0.6, 9 is the 
angle between the x and y coordinates: 9 = corresponds to the x axis of Figure [6] 




Figure 8. Dirichlet Region. The top and bottom faces in the z direction are identified 
with a rotation of 7r/3 



73 



74 



(\ 





a/2 


\ 




1 





yfla/2 









1 







J 


v 





1 


/ 




( V 2 




-VS/2 





o\ 


\/3/2 




1/2 
















1 


b 


V 













(22) 



which together with their inverses generate V. The translational part is generated by 
T = T/Zjq = {71,73,7!}, together with the inverses. The fundamental region, FR 
for T/Z 6 has volume V = 6^/3a 2 b and is given by 6 copies of the DR in Figure [HJ 
FR = {DR,-f A .DR,-flDR,-flDR,^iDR^lDR}. The identifications are given in 
Appendix A.3| 



With a reasoning analogous to the previous sections the following result is obtained 
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using ([7j) and 



+Z 
+Z 







9l 9y gl 





4 4 4 

91 91 91 








(4) v + Z 
(4) v + Z 



(4)*, + Z 








2* 4 





£y 



(4), 
(4), 



(4), 



(23) 



where 



t 



x 



+ y/V3 x-yjsfZ 



1/6 



-2x x — yy3 



,2/6 

a a 
-2x + 2y/VS -4y/V3 



a 

-x + yyS 



a 

x + y\^3 



3/6 



a 

2y/y/S 



4/6 



x 



+ y/Vz 



^fJbV 



a 

( ~ 2 

a 2 12 

\ 



,5/6 



a a 

2 



/2 \ 
a 2 



366 2 



(24) 



/ 



The values of a = 1.2, b = 1.4 are chosen. The truncation in is done such that 
the summation is over a geodesic sphere of radius R ~ 100 of integer numbers. The 
corresponding truncation error on the energy density should be less than 1%. The result 
is shown in Figure El and is independent of the z coordinate. 

3.6. * = Z 2 <g> Z 2 

The face pairing matrices are now 



7i 



73 
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a\ 
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-1 








o\ 
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-1 
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-1 
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-1 
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/ 













(25) 
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Figure 10. Dirichlet Region. The opposite faces in the x direction are identified with 
a rotation of tt. The adjacent faces are identified with a rotation of 7T, also. 



which together with their inverses generate T, and the volume of the manifold is V = 
2abc, Figure [TUl The translational part is generated by T = Y / (Z 2 ®Z 2 ) = {7^, 71, 7273}, 
with the inverses. The fundamental region, FR for T/(Z 2 <S> Z 2 ) has volume V = 8abc 
and is given by 4 copies of the DR in Figure TlO\ FR = {DR^i.DR^2-DR^i^.DR}. 
The identifications are given in Appendix A. 4 With ([7]) and ([9j) the Casimir energy 
results in 






9l 9l 9 2 z 






(4)„ + Z 



(4), 







9] gl 



(4), 





al gl gl 



(4), 



(26) 



where 
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-0.6 -0.4 -0.2 0.2 0.4 0.6 



X 

Figure 11. The Casimir energy density in natural units, given in (T3(5|) inside the DR. 
in Figure [TU1 a cylindrical region surrounding the x axis, of radius 0.762, is the angle. 



a — 2x b — 2y 



2a 

( 4a 2 


V 



,1/2 
\ 



(27) 



' 2b 



46 2 
4c 2 

1.4 and c = 1.3 are chosen. The truncation in is done 
such that the summation is over a geodesic sphere of radius R ~ 100 of integer numbers. 
The corresponding truncation error on the energy density should be less than 1%. The 
result is shown in Figure [TTJ This is the only manifold for which there is not an axis 
such that the energy density expectation value is constant. 



The values of a 



1.2, b 



4. Conclusions 



The Casimir energy for the oriented closed Euclidean space-times are obtained in this 
work. 

Closed manifolds are described by the discrete subgroups of the isometry group. 
The crystallographic groups are discrete subgroups of the full isometry group of 
Euclidean space. Among all crystallographic groups, only a few of them are of interest 
in connection to manifolds. These groups must be torsion-free, for the Euclidean and 
Hyperbolic spaces. In the celebrated list of problems proposed by Hilbert, the 18th was 
answered affirmatively by L. Bieberbach in 1910-1912 [23J, [24], see also [25]. Bieberbach 
showed that there are only a finite number of crystallographic groups for an Euclidean 
space of fixed dimension. 

When the dimension of the Euclidean space is 3 there are only 10 compact 
manifolds, of which 6 are orientable and 4 non orientable. 

In this present work we have investigated the 6 compact orientable manifolds and 
obtain the Casimir energy for each of them. Due to the more simple structure connected 
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to these topological spaces compared to the hyperbolic or spherical case, an analytical 
result is obtained in terms of Epstein zeta functions. It is explicitly obtained that 
our result is independent of the type of coupling with curvature £ for these particular 
manifolds, which means that the calculated Casimir energy density is the same even in 
an inflationary de Sitter regime. 

According to the theorem in Section [21 it is always possible to find the normal 
subgroup of finite index, T. With this, the n— torus that covers the manifold a finite 
number of times is obtained. The Casimir energy is given by a finite sum of Epstein 
type zeta functions. The number of distinct zeta functions is given by the index of T in 
T, which is the number of times the manifold is covered by this n— torus. 

A numerical plot inside the Dirichlet Region for each manifold is given. The 
manifold in Section 13.11 is the only one for which the Casimir energy density is 
homogeneous. For the manifolds Sections 13.21 - 13.51 the Casimir energy density depends 
on the x and y coordinates inside the fundamental domain. For the last manifold in 
Section 13.61 the Casimir energy density depends in all the coordinates. 

The sign of Casimir energy density is negative for all manifolds. Remind that there 
is a subtraction of the divergence of the infinite R 3 covering space, which means that 
the energy of the manifold is less than the energy of the infinite R 3 . The effect of 
choosing different representations for the fundamental group T acting on the field, like 
twisting, can increase the energy of the vacuum [22]. deWitt finds that the sign of the 
energy depends on the representation of the fundamental group, the orientability of the 
manifold and the spin of the field. This must have a more fundamental reason connected 
to the Atiyah-Singer theorem, and is beyond our knowledge. 

It is explicitly checked that the truncation error should be on the order of 1/R 
where R is the geodesic radius of the sphere of integer numbers where the truncation 
of (j7]) was performed. It was used R ~ 100 such that the truncation error should be on 
the order of 1%. The respective values of the parameters of the manifold are indicated 
in each case. 
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Appendix A. 

Appendix A.l. ^ = Z4 

r = (x, y, z) — > (x + ma, y + na, z + 4kb) 
— > (—y + ma, x + na, z + (4k + 1)6) 
— > (—x + ma, — y + na, z + (4k + 2)6) 
— ► (y + ma, —x + na, z + (4k + 3)6) 
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According to the above the Casimir energy is given by 

+ OC 

Pk = -~2 E ' {[{maf + {naf + {AkbfY 2 

m,n,k=— oo 

+ [(-x - y + ma) 2 + (x - y + na) 2 + ((4k + l)b) 2 ] ~ 2 

+ [(-2x + ma) 2 + (-2y + na) 2 + ((4k + 2)b) 2 ] ~ 2 

+ [(y-x + ma) 2 + (-x - y + na) 2 + ((4k + 3)b) 2 ] ~ 2 j 

Appendix A. 2. ^ = Z 3 

f — (x, y, z) — > [x + ma + na/2, y + naV3/2, z + 3kb] 
— > [—(x + yV3)/2 + ma + na/2, 
(iv^ -y)/2 + naV3/2, z + (3k + 1)6] 
— > [— (x — yV3)/2 + ma + na/2, 
- (xV3 + y)/2 + naV3/2, z + (3k + 2)6] 

According to the above the Casimir energy is given by 

+00 



°k = 2 E ' { [( ma + na / 2 ) 2 + ( n °v / 3/2) ; 

m,n,k=— oo 

+ (3A;6) 2 ] " 2 + [(-(3a: + j/v / 3)/2 + ma + na/2) 2 
+((V3x - 3y)/2 + nav 7 ^) 2 + ((3k + l)6) s 



+ 



-(3a; - yV3)/2 + ma + na/2f 



+ (-(V3x + 3y)/2 + naV3/2) 2 + ((3Jfe + 2)6) ; 



Appendix A. 3. * = Z 6 



f = (a:, y, 2;) — > [a: + ma + na/2, y + nay/3/2, z + 6/co] 

— > [(a; — yV3)/2 + ma + na/2, 
(v^x + y)/2 + na v^/2, z + (6fc + 1)6] 

-> [-(a; + yV3)/2 + ma + na/2, 

(iv^ -y)/2 + naV3/2, z + (6k + 2)6] 

[— x + ma + na/2, — y + naV3/2, z + (Qk + 3)6] 

— > [— (x — y\/3)/2 + ma + na/2, 

- (xv^ + y)/2 + na >/3/2, z + (6k + 4)6] 
— > [(x + yV3)/2 + ma + na/2, 

- (v^x - y)/2 + naV3/2, z + (6k + 5)6] 
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According to the above the Casimir energy is given by 

j +00 

Pk = -^2 E '{[(™ + m/2) 2 + M/2) 

m,n,k=— 00 

+ (6fc6) 2 ] " 2 + [(-(x/2 + yV3)/2 + ma + na/2) 2 



+((V3a; - y)/2 + naV3/2) 2 + ((6k + l)bf 



-2 



+ 



(-(3a; + yVs)/2 + ma + na/2) 



+((v / 3x - 3j/)/2 + nav^^) 2 + ((6fc + 2)6)' 



+ 



-2x + ma + na/2) 2 + (-2y + naV3/2) 2 
+ ((6k + 3)6) 2 ] " 2 + [(-(3a; - yv / 3)/2 + ma + na/2f 
+ (-(V3x + 3y)/2 + naV3/2) 2 + ((6Jfe + 4)6) 2 



+ 



;-(x - y>/3)/2 + ma + na/2f 



+(-(V3a; + y)/2 + naV3/2) 2 + ((6k + 5)bf 

Appendix A. 4- ^ = Z 2 ® Z 2 

r = (x, 2) — > (x + 2ma, y + 2nb, z + 2kc) 
— > (x + a + 2ma, — 7/ + 2no, —2; + 2/cc) 
— ► (—x + 2ma, ?/ + b + 2no, — z + c + 2/cc) 
— > (—x + a + 2ma, — y + 6 + 2no, z + c + 2/cc) 

According to the above the Casimir energy is given by 

+00 



~ '{ P ma ? + ( 2nb ) 2 + ( 2kc ) 2 ] 



m,n,k=—oo 



+ [(a + 2ma) 2 + (-2y + 2nb) 2 + (-2z - 2/cc) 2 ] 2 
+ [(-2a; + 2ma) 2 + (b + 2n6) 2 + (-2z + c + 2/cc) 2 ] ~ 2 
+ [(-2x + a + 2ma) 2 + (-2y + 6 + 2n6) 2 + (c + 2/cc) 2 ] ~ 2 
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